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Abstract 

Stress measurement methods using X-ray diffraction (XRD methods) are based on so-called fundamental 
equations. The fundamental equation is described in the coordinate system that best suites the measurement 
situation, and, thus, making a comparison between different XRD methods is not straightforward, ffow- 
ever, by using the diffraction vector representation, the fundamental equations of different methods become 
identical. Furthermore, the differences between the various XRD methods are in the choice of diffraction 
vectors and the way of calculating the stress from the measured data. The stress calculation methods can 
also be unified using the general least-squares method, which is a common least-squares method of multi¬ 
variate analysis. Thus, the only difference between these methods turns out to be in the choice of the set 
of diffraction vectors. In light of these ideas, we compare three commonly used XRD methods: the sin 2 1 // 
method, the XRD 2 method, and the cos a method using the estimation of the measurement errors. 
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I. INTRODUCTION 


The cos a method [ l|], an X-ray diffraction (XRD) method, is widely used in industry, but there 
are few studies comparing it with other XRD methods in their theoretical aspects. Although we 
gave a mathematical explanation of the methods based on Fourier series for the plane stress (biaxial 
stress) case [2], it is important to place the cos a method and the other XRD methods on a common 
mathematical basis. In this study, we compare the cos a method for the triaxial stress case [3y with 
the sin" ip method (for example, please see [4j]) and the XRD" method [5] from the aspect of the 
fundamental equation. First, we show that all three methods are based on a common fundamental 
equation in the diffraction vector representation. Second, we show that this fundamental equation 
can be solved in a common way by using the general least squares method [ 6 ]. Accordingly, the 
only difference between XRD methods is the choice of the set of diffraction vectors. Finally, we 
compare XRD methods based on the measurement error estimation. 


II. FUNDAMENTAL EQUATION 


For the sake of simplicity, we will suppose that the specimen is a polycrystal composed of 
elastically isotropic crystallites. Furthermore, we will assume that the microscopic stress of the 
specimen can be ignored. 

We will use the conventional coordinate system of the sin 2 ip method (for example, see Fig. 2 
of [4]). The unit diffraction vector n (in the following, we call the unit diffraction vector the 
“diffraction vector”) can be described by two angles: (p and ip (Fig. 1). (p is the rotation angle of 
the diffraction vector around the S 3 axis, and ip describes the tilt angle of the diffraction vector 
from the S 3 axis. Though [4] describes the strain corresponding to this diffraction vector as s 
we will consider diffraction by a single diffraction plane (hkl) and use for simplicity. The 
X-ray measured strain can be described using the strain in the specimen frame of reference as 


= g| 1 cos 2 (p sin' \p + £22 sin 2 (p sin' ip + £ 33 cos 2 ip 
+ £12 sin( 20 ) sin 2 ip + s l3 cos (p sin( 2</9 
+ £23 sin cp sin (2 ip) 


( 1 ) 


This is the fundamental equation of the sin' ip method (for example, Eq. (13) of [4]). Accordingly, 


the sin 2 ip method can be considered an inverse problem of estimating £ (/ (z, j = 1 • • • 3) from s, 
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measured with a certain set of (< f, if/). 



FIG. 1. Arrangement of X-ray stress measurement. 


The diffraction vector n can be described using {(p, ip): 



f \ 

n\ 


( \ 
sin ip cos (p 

n = 

n 2 

= 

sin ip sin (p 




, cos Ip J 


Substituting Eq. ([2]), Eq. ([[]) becomes 

e n = n\e n + n\s 2 2 + n\s 33 + 2n i n 2 sn + 2n l n 3 s l3 + 2n 2 n 3 s 23 


( 2 ) 


( 3 ) 


This is the fundamental equation in the diffraction vector representation. 

As in the case of Eq. ([2]), a diffraction vector n can be represented by two circumference 
angles. Here, we will represent n by a (cp, ip) pair, which is equivalent to the (cp, ip) pair of the 
sin 2 ip method (Fig. 1). Using (0, ip) pairs, the diffraction vector can be displayed in a pole figure. 
Figure 2a shows the definition of the pole figure (angles are in radians). This figure displays a 
diffraction vector as a (0, ip) pair and shows the set of diffraction vectors as a constellation. The 
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center of the figure is (<p = 0, ip = 0), and d> is the circumference angle. The distance from the 
center describes sin < p. 

Figure 2b shows an example of a diffraction vector of the sin ip method (cf) = 0, tp = 45°). It 
has to be emphasized that when measuring s n for a tp) pair with a position-insensitive X-ray 
detector, which many instruments of the sin 2 tp method use, several irradiations and detections 
are required in order to find the peak position of the diffraction ring. On the other hand, XRD 
instruments with an area detector require only one X-ray irradiation and detection to find the peak 
position. 

(a) (b) 


> = 0 


(]) = 71/2 


y = 0 

y = n/6 
y = 7i/3 
y = 7i/2 


= 3ti/2 


FIG. 2. (a) Definition of the pole figure (angles are shown in radians), (b) Example of a diffraction vector 
of an X-ray irradiation of the sin 2 <p method corresponding to Up - 0, tp - 45°). 


A. cos a method 


The cos a method measures the stress from one or more Debye-Scherrer (D-S) rings. Figure 2 
in [lfl, Fig. 1 in [7] or Fig. 1 in [2] illustrate the set up of this method. The diffraction vector of the 
cos a method ( Eq. (5) of [1]) is 



( \ 
n i 


( \ 
cos T] sin ipo cos 0 O - sin rj cos <po cos 0 O cos a - sin ij sin tpo sin a 

n = 

n 2 

= 

cos T] sin \pQ sin (p 0 - sin tj cos ipo sin 0 O cos a + sin 77 cos (p 0 sin a 




cos tj cos ipo + sin r/ sin ip 0 cos a 


( 4 ) 
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Note that ip 0 of this method is not identical to \p of Fig. 1. Using Eq. ©, the fundamental equation 
of the cos a method (for example, Eq. (8) of [!]) becomes identical to Eq. ©. 

The diffraction vector of Eq. © can be expressed as an equivalent (0, ip) pair of the sin 2 \p 
method, as follows: 


ip = cos 1 n 3 


(5) 


and 


<P = 


tan 1 (n 2 /ni) 


(if n 2 > 0) 


( 6 ) 


tan l (n 2 /n l ) + n (if n 2 < 0) 


Figure 3a shows an example pole figure of a constellation of diffraction vectors resulting from an 
X-ray irradiation. The conditions of the figure are taken from [2]: 26 = 157.08° (r/ = 11.46°), 
0o = 0°, ipo = 35°, and 0° < a < 360°. It has to be emphasized that this constellation corresponds 
to a single X-ray irradiation. Because the cos a method utilizes the data from a whole D-S ring, 
it is possible to measure the biaxial stress with a single X-ray irradiation. To measure the triaxial 
stress, the cos a method requires a number of X-ray irradiations with two to four (0 O , ip 0 ) pairs 

d so. 


(a) 


(b) 




FIG. 3. (a) Example pole figure of the diffraction vector of the cos a method for an X-ray irradiation, (b) 
Example pole figure of the diffraction vector of the XRD 2 method For an X-ray irradiation. 
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B. XRD 2 method 


The XRD 2 method [5] measures the stress from fractions of D-S rings. To discuss this method, 
we will use the coordinate system depicted in Figs. 6 and 10 of [9]. The diffraction vector n = 
(n u n 2 , n 3 ) T of the method (Eq. (5) in [9] and shown as (hi, h 2 , h 3 )) is 


n i = sin @(sin tp sin tp sin to + cos 0 cos to) + cos 8 cos y sin tp cos tp 
- cos 8 sin y(sin tp sin tp cos to - cos tp sin to) 
n 2 = - sin 0(cos <p sin tp sin to - sin tp cos to) - cos 8 cos y cos tp cos tp 
+ cos 8 sin y(cos tp sin tp cos to + sin tp sin to) 
n 3 = sin 8 cos tp sin to - cos 8 sin y cos tp cos to- cos 8 cosy sin tp (7) 


The fundamental equation of the XRD 2 method in the diffraction vector representation is identical 
to Eq. ©. 

The diffraction vector of Eq. ([7]) can be expressed as an equivalent (tp, tp) pair of the sin 2 tp 
method using Eqs. © and ©. Figure 3b shows an example of a constellation of diffraction 
vectors of the XRD 2 method resulting from an X-ray irradiation. To make the difference from the 
cos a method clear, the angles are: 28 = 157.08°, tpo = 90°, tp = 35° (this tp is not identical to that 
of Fig. 1), to = 90°, and 62.5° < y < 117.5°. The range of y was taken from [K)|]. Comparing 
Figs. 3a and 3b, it can be seen that the constellation of diffraction vectors of the XRD 2 is part of 
that of the cos a method. Thus, the cos a can measure the stress by using less X-ray radiation than 
that of the XRD 2 method. 


C. Comparisons of diffraction vector formulas 

So far, we have seen that the fundamental equations of the sin 2 tp method, the cos a method, 
and the XRD 2 method are identical in the diffraction vector representation. In this section, we 
demonstrate that the expressions of the diffraction vectors (i.e. Eqs. ©, ©, and ©) agree each 
other with a proper coordinate transformations. First, we show that Eq. © is a special case of Eq. 
©. Then we show that Eq. © is a special case of Eq. ©. 

sin 2 tp can be regarded as a method that measures only one point on a D-S ring: a = 0 of the 
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cos a method. Thus, substituting n,\ of Eq. © with a = 0 and cp 0 = cp, we obtain 


77 1 = cos rj sin ip 0 cos cp 0 - sin 77 cos 1 J/q cos cpo cos a - sin 77 sin cp 0 sin a 


= cos 0(sin ijjQ cos tj - cos ip 0 sin 77) 
= cos cp sin((^ 0 - 77 ) 


Using ipo = 1/7 + 77 (Fig- 2 of [lj), we obtain 


77 1 = cos cp sin ip 

With the same substitutions: a = 0, 0 O - 0, and t^o = (A + we find that Eq. © is equivalent to 
Eq. ©. Thus, the representation of the diffraction vector in the sin 2 1/7 method is a special case of 
that of the cos a method. 

Comparing the arrangement of the XRD 2 method with the arrangement of the cos a method, 
we find that 7 = n - a and 77 satisfies 9 = n/2 - 77. Thus, 771 of Eq. © can be modified as 


771 = cos ? 7 (sin cj) sin 1/7 sin u> + cos 0 cos m) - sin rj cos a sin 0 cos 1/7 
- sin tj sin cKsin (j> sin 1/7 cos c 0 - cos (p sin co) 


Furthermore, by setting co = n/2, cp = <p 0 + tt/ 2, and if/ = <// 0 , we obtain 


77 1 = cos 77 sin if /0 cos 0 O - sin 77 cos i/7 0 cos 0 O cos a - sin 77 sin sin a 

which is identical to 77 1 of Eq. ©. In the similar manner, Eq. © becomes identical to Eq. @ with 
the conversions: 

7 = n - a 
6 = n /2 - tj 
' CO - 7r/2 
(p = (po + tt /2 

ip = ip q 

Thus, the representation of the diffraction vector of the cos a method is a special case of that of 
the XRD 2 method. 


III. GENERALIZED STRESS DETERMINATION 

XRD methods can be regarded as inverse problems to obtain the strain of the specimen as 
the coefficients of Eq. © for a certain set of diffraction vectors. In the strict sense, the sin 2 ip 
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method and the cos a method solve the fundamental equation by using simplified analyses that Or- 
tner named “linear-regression methods” [11]]. Though linear-regression methods are useful when 
computational power is limited, they are not proper least-squares methods. The generalized least- 
squares methods of multivariate analysis, which directly solve Eq. ©, have been discussed by 
[6, Jil 12]. To make a simple comparison of the methods, we solely use the generalized analysis 


in the following. As [130 called the sirr i// method with the general least-squares method analy¬ 


sis the “generalized sin 2 1 \jj method”, we will call the cos a method with the general least-squares 
method analysis the “generalized cos a method”. We will not discuss the difference between the 
linear-regression methods and the generalized least-squares methods any further. 

Let us consider the case of observing s n with a set of k diffraction vectors. n ; = (n,i, n i2 , n i3 j' 
describes the z-th diffraction vector, and corresponding equivalent (<A, i/z,) pairs can be calculated 
using Eqs. © and ©. From Eq. ©, s n for n, satisfies 


= n 2 n e ii + Z7 2 2 e 22 + iij 3 s 3 3 + 2n n n i2 s l2 + 2n n n i3 s ]3 + 2n a n i3 s 23 


( 8 ) 


Let us define a k x 6 matrix embodying the coefficients of Eq. ® : 


F = 


77 1 1 77 12 77 13 2/7h/7i 2 2/7 n/713 2/7, 2 /7, 3 
77 2 , « 22 « 23 2n 2l n 22 2/7 2 j/7 23 2n 22 n 23 


\ n h n h n k 3 2n ki n k .2 2n H n k3 2n k2 n k3 


(9) 


Using this matrix, the set of fundamental equations can be described as 


t \ 

^ni 


"ri 2 


\ £ n k ) 


= F 


( \ 
£n 

^22 

^33 

^12 

£13 

V £ 23 j 


( 10 ) 


The strains (eu, s 22 , s 33 , s\ 2 , e ]3 , e 23 ) T can be related to the stresses cr = (crn, cr 22 , cr 33 , cr 12 , cr 13 , cr 23 )' 1 
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as 




1 -V -V 0 0 0 


cr a 

£22 


-v 1 -v 0 0 0 


CT 22 

s 33 

1 

-v -v 1 0 0 0 


c 33 

£12 

“ E 

0 0 0 1 0 0 


O' 12 

£13 


0 0 0 0 1 0 


O' 13 

£23 j 


y 0 0 0 0 0 1; 


,0-23; 


where E and v are the X-ray Young’s modulus and Poisson’s ratio, 
(fill) . Eq. (flQl) becomes 


Srr 


( 11 ) 


respectively. Substituting Eq. 


e n - F ■ Srr = Mcr 


( 12 ) 


The general least-squares solution of Eq. (fT2l) is 


rr - M^£ n 


(13) 


where M~ 1 ' is the Moore-Penrose’s general inverse of M. Equation (fl3l) is the universal solution 
for the XRD measurement. 

Using / n)j, the (/, j )th component of M\ and s nj , each component of <r can be described as a 
linear combination. For example, 

k 


O-U = 'Yj m \i s n i 
7=1 

If the measurement error of each s ni is independent and has a deviation 6s, the measurement error 
of cr n is 


dcrii 



(14) 


The errors of the other components of er can be estimated in a similar way. However, not all of 
the measurement points of the cos a method are independent, and the assumption that errors are 
independent is not fully satisfied. In this case, Eq. (fl4l) underestimates the error and an adequately 
sparse set of s ni is required to estimate the correct error. 


IV. COMPARISON OF TRIAXIAL STRESS MEASUREMENT 

In the previous sections, we described the way to calculate stress from s n measured for a set 
of diffraction vectors n. Though the stress can be calculated using Eq. (fl3l) . the equation does 
not tell which set of diffraction vectors should be chosen to measure the stress. However, once the 
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set of diffraction vectors is chosen, we can estimate the error of the stress measurement by using 
Eq. (fl4l) . In this section, we compare XRD methods in terms of their errors as estimated by Eq. 
(fl4l) . Specifically, we compared representative constellations of the sin 2 p method and the XRD 2 
method, three constellations in [3] and a new constellation of the cos a method. We assumed that 
the hkl = 211 diffraction plane of an cr-Fe specimen was measured with Cr-K a characteristic X- 
rays. The diffraction angle was taken to be 6 = 78° (i.e. rj = 12°), and X-ray Young’s modulus and 
Poisson’s ratio were 


E = 221 (GPa) 
v = 0.28 


Table 1 shows the (cp, p) pairs for the generalized sin 2 p method [14]. This set requires 31 


(p,p) pairs. In the case of the sin 2 p method, this means 31 individual data acquisitions (for 
convenience, we call them “frames” hereafter) are required. As stated before, when using a stress 
measurement instrument with a position-insensitive X-ray detector, one frame requires several 
X-ray irradiations. The pole figure of this constellation is shown in Fig. 4a. 


TABLE I. Constellation of {(p, p) pairs for the triaxial stress measurement with the sin 2 p method. 


<p 

p 

o 

O 

0°, ±18°, ±26°, ±33°, ±39°, ±45° 

45° 

±18°, ±26°, ±33°, ±39°, ±45° 

80 

O 

o 

±18°, ±26°, ±33°, ±39°, ±45° 


Table 2 shows the pairs of (< p , p) pairs for the XRD“ method. Note that the pairs ( (p , p) of this 
table are those of Eq. ([7]) and are not identical to the equivalent (cp, p) pairs of the sin 2 p method. In 
the following calculations, we set <x> - 110° and 70° < y < 110° (in 5° step) [15Q. This set consists 
of 33 frames (data acquisitions). In the case of the XRD 2 method, one frame can be acquired with 
a single X-ray irradiation. The pole figure of this constellation is shown in Fig. 4b. 

Table 3 shows the (0 O , Po) pairs for the generalized cos a method. Type A is according to 
[8], and Types B and C are according to [3]. Type D is new. Type A requires two frames (data 
acquisitions), Type B requires four frames, and Types C and D require three frames. Compared 
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FIG. 4. (a) Example pole figure of the diffraction vector of the sin 2 1 // method for a triaxial stress measure¬ 
ment. (b) Example pole figure of the diffraction vector of the XRD 2 method for a triaxial stress measure¬ 
ment. 


TABLE II. Constellation of (<f>, ifi) pairs for the triaxial stress measurement with the XRD 2 method. 


<p 

<A 

o 

O 

o 

O 

15° 

0°, 45°, 90°, 135°, 180°, 225°, 270°, 315° 

30° 

0°, 45°, 90°, 135°, 180°, 225°, 270°, 315° 

45° 

0°, 45°, 90°, 135°, 180°, 225°, 270°, 315° 

O 

O 

0°, 45°, 90°, 135°, 180°, 225°, 270°, 315° 


with the other methods, the cos a method requires the fewest frames. Figures 5a-d show the pole 
figures of Type A-D. 

The generalized cos a method calculates the stress using whole D-S rings. The number of data 
points of one frame is n a = 500 [ 2|]. On the other hand, the error of the stress her estimated using 
Eq. (fl4l) is proportional to 1 / yfn^. From this, one may conclude that the accuracy of the stress 
measurement can be infinitely improved if n a is increased. But as stated previously, the neigh¬ 
boring points of a frame are correlated with each other and the effective number of independent 
data points is less than 500. Here, we will not discuss the most proper n a , but will instead assume 
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(a) 


(b) 


o 

o 





o 

r\ r~\ 

U U 



FIG. 5. (a) Pole figure of the diffraction vector of the cos a method for a triaxial stress measurement (Type 
A), (b) Pole figure of the diffraction vector of the cos a method for a triaxial stress measurement (Type B). 
(c) Pole figure of the diffraction vector of the cos a method for a triaxial stress measurement (Type C). (d) 
Pole figure of the diffraction vector of the cos a method for a triaxial stress measurement (Type D). 


n a = 72 (5° step) in accordance with the XRD 2 method. This assumption is realistic for the error 
estimation and sufficient for the purpose of comparison with other methods. 


Table 4 shows the error estimated using Eq. (fT4l) . Though the values of the sin 2 tfj method are 
not identical to those of [6], the differences are small. The reason for these small discrepancies 
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TABLE III. Constellations of (<f>o,tj/o) pairs for the triaxial stress measurement with the generalized cos a 
method. 


Method name 

Combinations of (0o, 1 A 0 ) 

No. of Frames 

Type A 

(0°, 45°), (180°, 45°) 

2 

Type B 

(0°, 45°), (90°, 45°), (180°, 45°), (270°, 45°) 

4 

Type C 

(0°, 45°), (90°, 45°), (0°,0°) 

3 

Type D 

(0°, 45°), (120°, 45°), (240°, 45°) 

3 


is under investigation. The XRD 2 method which consists of 33 frames showed the best accuracy. 
Compared with the generalized sin 2 if/ method, the XRD 2 method is approximately six times more 
accurate and uses a similar number of frames (31 frames). 


TABLE IV. Comparisons of the error estimated using Eq. ( 1 1 41 ). be = 10 4 and n a = 72 were assumed. 



Estimated errors (MPa) 


Scr ii 

her 22 

Sc 33 

Scr 12 

Scr 13 

Sc 23 

sin 2 (7 method 

36.6 

36.6 

15.4 

26.0 

7.0 

7.0 

XRD 2 method 

5.6 

5.6 

2.8 

3.6 

2.0 

2.0 

cos a method 

Type A 

50.7 

138.1 

50.7 

8.4 

1.8 

8.4 

Type B 

6.2 

6.2 

3.1 

5.9 

1.8 

1.8 

Type C 

15.7 

15.7 

8.3 

9.4 

5.8 

5.8 

Type D 

8.2 

8.2 

3.6 

6.0 

3.3 

3.3 


The generalized cos a method showed good accuracy when more than three frames are taken 
(i.e., Types B-D). Type B with four frames is as accurate as the XRD 2 method. This result can 
be understood intuitively in that a single frame of the XRD 2 method takes l/8th of the D-S ring, 
while a single frame of the cos a method acquires a whole D-S ring. Thus, the cos a method 
can achieve similar accuracy with l/8th of the frames of the XRD 2 method. Moreover, by using 
Type D, we can reduce the number of frames by one while losing only lose 30% of the accuracy. 
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Consequently, we recommend Type D for the triaxial stress measurement with the generalized 
cos a method. 

V. SUMMARY 

This study showed that the sin 1 2 3 4 5 6 7 ijj, cos a, and XRD 2 methods can be described with a common 
fundamental equation using the diffraction vector representation. By fitting the data with the gen¬ 
eralized least-squares method, the only differences between these methods are in the choice of 
the set of diffraction vectors. The differences between the sets of diffraction vectors become clear 
in the pole figure plot. We also estimated the errors of the XRD methods for typical choices of 
diffraction vector and demonstrated that the XRD 2 method with 33 frames is the most accurate. We 
further showed that the generalized cos a method with four frames is comparable in accuracy to 
the XRD 2 method. However, from the viewpoint of the balance between the number of the frames 
and the accuracy, the generalized cos a method with three equally spaced frames is recommended. 
In the future, the authors will test the conclusions of this study by making actual measurements. 
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